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Abstract 



We construct the gravity background which describes the dual field theory with 
aging invariance. We choose the decay modes of the bulk scalar field in the inter- 
nal spectator direction to obtain the dissipative behavior of the boundary correlation 
functions of the dual scalar fields. In particular, the two-time correlation function 
at zero temperature has the characteristic features of the aging system: power law 
decay, broken time translation and dynamical scaling. We also construct the black 
hole backgrounds with asymptotic aging invariance. We extensively study characteris- 
tic behavior of the finite temperature two-point correlation function via analytic and 
I/-) numerical methods. 
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1 Introduction 

AdS/CFT correspondence has been successful in describing strongly-coupled field theories 
using weakly-coupled classical gravity backgrounds, providing new paradigm of analytical 
methods for theories with strong coupling pQ [2] . Recently this correspondence extended its 
application to non-relativistic setup with general dynamical exponent z ^ 1, Those non- 
relativistic holographic theories include Schrodinger [21IHG3I6] an d Lifshitz holography [7]j^] 
Since then, they have been enjoying continuous attractions. 

An acquaintance of non-relativistic holographic theories with z ^ 1 tells that we can 
construct a plethora of explicit time dependent backgrounds by utilizing the distinguished 
role of time. This is different from the relativistic theories. Thus we might hope to shed 
lights on time dependent backgrounds in string theory and in AdS / CFT correspondence even 
though they are in general known to be challenging. We would like to add some contributions 
to this direction in the context of the non-relativistic Schrodinger holography with z=2. 

One of the simplest applications of the non-equilibrium physics is known as aging. Aging 
phenomena can be characterized by two-time correlation functions, which typically show 
that older systems relax in a slower manner than younger systems after entering a given 
physical phase under study, see e.g. [HJPH] for reviews. Two pioneering works of aging in 
holographic setup have been put forward in [TT] [12] by generalizing Schrodinger background 
[3 J |4J with explicit time dependent terms. The authors of p2] observe that aging algebra can 
be realized as a subset of Schrodinger algebra using a singular time dependent coordinate 
transformation. See also [13] for similar observations of general time dependent deformations 
of Schrodinger backgrounds. 

Along the line of Schrodinger backgrounds [2] [I], there have been progresses with much 
interest. Their finite temperature generalizations with the thermodynamic analysis can be 
found in [HI [151 EH HZj and transport properties are analyzed in [18] , recently. These so- 
lutions are constructed nicely using the null Melvin twist [19|[20j. There also have been 
efforts constructing Schrodinger solutions from string or M theory, sometimes with super- 
symmetry. See [211 E2J E31 EH EH] , for example. While these Schrodinger backgrounds have 
been drawn much interest, the progresses directly along this line have been hindered due 
to the conceptual difficulties related to their non-trivial boundary structures. Furthermore 

^ee e.g. [5] for a very different usage of the Lifshitz symmetry with z=3 in the context of quantum field 
theory of gravity. 
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one quickly finds that finite temperature generalizations of the Schrodinger backgrounds are 
complicated, and thus practical calculations are difficult. 

An alternative candidate for the Schrodinger holography, AdS in light-cone, has been 
constructed in [5][6J even before the finite temperature generalizations of the Schrodinger 
backgrounds. Its finite temperature generalizations were done in [I5][26] with thermody- 
namic analysis, and transport properties are also analyzed in [26] [27]. It turns out that 
these two different holographic theories, Schrodinger backgrounds and AdS in light-cone, 
have the same thermodynamic properties [H][2S] and also the same transport properties 
when the comparisons are reliable [26]. Thus AdS in light-cone is a viable candidate for 
Schrodinger holography. While it remains as a question whether these two theories describe 
the same physical properties or not, AdS in light-cone has advantages over the Schrodinger 
backgrounds: it has well-defined holographic renormalization and is as simple as original 
AdS in its computations. 

Motivated by these observations, we would like to construct a holographic model of 
aging in the context of [5] [6], Aging in light-cone, by realizing the aging symmetry with 
the singular time dependent coordinate transformation. This is done in section [2] From 
the bulk gravity side, we obtain the two-point correlation function of the dual field theory. 
In particular, we consider the, so-called, spectator coordinate internal direction]^] 

Furthermore, since we are dealing with time dependent and dissipative system, we choose the 
decay modes of the bulk scalar field in the internal spectator direction. This is in contrast 
with the complexification of spacetime done in [T2]. As a result, we obtain generalized two- 
time correlation function of the dual field theory, which exhibit the characteristic features of 
the aging system. 

In section 3 we consider finite temperature generalizations of the aging holography on 
Aging in light-cone as well as Aging in Schrodinger backgrounds. We focus on the former 
case, which has simpler form, and then give some comments on the latter. We study the 
finite temperature two-point correlation function in Aging in light-cone via AdS/CFT cor- 
respondence. We provide the novel time dependent structure of aging two-point correlation 

2 Note that the geometric realizations of Schrodinger symmetry separates itself from the other holographic 
examples, providing a correspondence between (d, l)-dimensional field theory and (d+3)-dimcnsional gravity 
theory. In addition to the usual radial direction, it has another coordinate x~, which does not contribute 
to the coordinates in field theory. The x~ coordinate is known to provide a non-relativistic particle number 
or mass, referred here by M. We call this spectator direction, x~ , as internal because it does not show up 
explicitly in the dual field theory. 

While the full potential associated with the x~ direction remains to be seen, we believe that the treatment 
adapted in this paper is one of the first non-trivial examples of utilizing this spectator direction. See e.g. 
28J[29 for different attempts on this issue. 
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function in coordinate space, compared to that of Schrodinger case and to that of momentum 
space. As the latter is not exactly solvable, we use numerical as well as analytic methods 
to see its behavior. Finally, we give some comments on the aging phenomena from the 
Schrodinger background stressing the difference compared to those of the Aging in light- 
cone. It turns out that the internal isometry direction plays a crucial role throughout in this 
paper. In section 4, we give conclusion. 

2 Zero Temperature Aging Holography 

In this section, we start with a brief review of aging holography based on recent papers 
[TT][T2] and also fix some notations. The relevant geometry is simple and possesses lots 
of isometries, yet there are many difficulties in its fundamental level. We also notice [T3] , 
which considered some universal time dependent deformations in the context of Schrodinger 
geometry. One of the deformations is essentially the same as that of [12], yet the discussion 
was brief and not developed to the degree done in [T2"] . 

Then we develop the basic two-point correlation functions from a background called AdS 
in light-cone [HI EJ [15J, [26] following [TT] [T2] . On the way, we develop several conceptually 
important points of the applications of AdS/CFT on time dependent backgrounds. They 
include 

i) Establishing the holographic aging equation which describes non-equilibrium critical phe- 
nomena, 

ii) Non-trivial radial fall offs of wave function, which are different from the conformal dimen- 
sions of the corresponding field, and 

iii) Utilizing the x~ coordinate as an internal direction from the point of view of holography, 
which gives an alternative way to get real two-time correlators instead of complexifying the 
geometry done in [T2] . 

To make a clear distinction between these two approaches, we refer Aging background 
and Schrodinger background for the geometric realization with the additional terms in the 
background [3] [I] [12], while we use the term Aging in light-cone and AdS in light-cone for 
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2.1 Zero Temperature Aging Background : A Review 



Starting with the Schrodinger background [3] [I] 



ds 2 = r 2 (dy 2 - 2dx + dx~ - /3 2 r 2 dx +2 ) + , (1) 



the following Aging background is constructed in [12 



ds 2 = r 2 



df - 2dx + dx- - ( P 2 r 2 + °^ ) dx +2 



dv 2 

2a/3rdrdx + + — . (2) 



This turns out to be one of the simplest time dependent background which possesses a large 
number of isometries called aging symmetry with dynamical exponent z=2. The Dilatation 
scaling symmetry D is associated with the transformations, y — > Ay, x + — > X 2 x + and x~ — > 
x~ . Other relevant symmetries are spatial translations Pi, Galilean boosts Ki, rotations Rij, 
central element M interpreted as a particle number and the special conformal transformation 
C , while the time translational symmetry is broken by the explicit time dependence of the 
metric globally. We come back to this below. 



There exist several valuable observations in [T2], which are related to the geometry in 
equation ^ and worthwhile to be mentioned: 

• The metric ^ is connected to the Schrodinger background Q by a local, but singular, 
coordinate transformation 

x~ — > x~ — In (r 2 x + ) . (3) 

This singular time dependent coordinate change is a key to realize aging symmetry. 
It does not generate curvature singularities at any point in spacetime. It is argued in 
[T2] that this coordinate singularity in this context is similar to a black hole horizon, 
where physical boundary conditions are imposed, and thus has profound effects, from 
the gauge/gravity point of view. This geometry is called locally Schrodinger in [T2]. 
Thus it is important to impose physical boundary conditions on the time boundaries 
in addition to the spatial boundaries. 

• Thus, locally, the Aging background has the full Schrodinger symmetry with modified 
set of generators, especially the time translation (Ha) and special conformal transfor- 
mation (Ca) generators as 

Ha = d x + — -z^-d x - , 
2x + 

h2o , 1 f-fi 1 



C A = -x + rd r + x + y-d + x +2 d x+ + - I £ * + ^ - af3x + ) d x - . (4) 
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The other relevant symmetry generators are not changed by the coordinate trans- 
formation ^ [12]. Globally the time translation symmetry is broken and the aging 
symmetry is realized as conformal Schrodinger symmetry modulo time translation, and 
described by the generators {Rij, Pt,Ki,D,C,M}. 

• The analysis is done with the complex x + , especially at x + = 0, and also the geometry 
is extended to be complex with complexified parameter a in [12]. With real a, the 
resulting correlator has the time dependence only in its phase. To achieve a relaxation 
process, a should be complex. 

We show that there is a way out from complexifying the geometry by utilizing the 
fact that x~ can be considered as an internal coordinate from the holographic point of 
view. 



2.2 Aging in Light-Cone 

In a similar manner, we start with the AdS in light-cone [U E] 

dr 2 

ds 2 = r 2 (dy 2 -2dx + dx-) + — , (5) 

which is shown to be another viable candidate for the Schrodinger holography. This metric 
can be obtained from the AdS metric with the following coordinate change 

1 

x + = b(t + x) , x~ = — (t — x) . (6) 

This coordinate transformation was introduced in [15j[26j[27] and should be viewed as a 
two-step procedure: a boost in the x-direction with rapidity log b, followed by transforming 
to light-cone coordinates. To ensure z=2, we assign [b] (the scaling dimension of b in the 
unit of mass) as —1, and thus [yi] = — 1, [x + ] = — 2 and [x~] =0. To achieve the desired 
Schrodinger isometry, it is further required to identify the coordinate x + as time and also 
to have a momentum projection along the other coordinate x~ . x~ direction provides an 
isometry to the background and thus the corresponding momentum gives the central element 
M which serves as the total particle number or the mass of the resultant non-relativistic 
theory. Furthermore this AdS in light-cone is known to have a well defined holographic 
renormalization and as simple as the case of AdS black holes [26] |^] 

It has been known that the Schrodinger background and the AdS in light-cone share the 
same physical properties such as thermodynamic [H)J [25] and transport properties [26J . We 

3 See e.g. |30j |31j for holographic renormalization in general, [32] for Schrodinger background using mod- 
ified definition of the stress tensor, and |33j for general anisotropic backgrounds using anisotropic scaling. 
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would like to investigate whether this similarity extends to this time dependent settings or 
not. 



With the singular coordinate transformation (|3j) and a notation change a(3 — > a, we 
obtain the following metric 



/ ol \ dr 2 

ds 2 = r 2 ( dy 2 — 2dx + dx~ dx +2 ) — 2a rdrdx + H — , (7) 

V x + J r l 

which is simpler than that of the Aging background ([2]). One may note that a is a dimen- 
sionless parameter. It can be shown that this metric shares the similar properties as the 
geometry ([2]), mentioned in the previous section 2.1 First, the metric ([7]) is locally AdS in 
light-cone, while time translation symmetry is broken globally by the singular coordinate 
transformation which connects these two metrics ([7]) and (J5|. Second, the isometry of the 
metric ([7]) is aging symmetry with the same set of the generators given in Q. 

With these properties in mind, let us investigate the background For this purpose, 
we change the coordinate as u = — and concentrate on the five dimensional background 
with dy 2 = dy\ + dy\ and the metric becomes 

ds 2 = \dy\ + dy\ — 2dx + dx~ — dx +2 H dudx + + du 2 

u 2 \ x + u 

where the boundary sits at u = 0. Let us couple the background with a probe scalar 0, 
whose action has the following form 

S = K [d 4 x [ duj—g (g uu dj> d u <\> + g ux " dj> d x -<\> 

J J Ug 

+g ux ~ d x -4> d u <j) + <T d v <f> + m 2 0) , (9) 

where K = — 7t 3 L 5 /4k 2 and /i, v = +, — , y, z. We used ub for the boundary cutoff, which is 
small and will be sent to zero eventually. The linearized field equation for becomes 

2M (i 9 + aM \ 
dx + 2x + J 

d 2 6 , Add) [AiMa + a 2 . 

Z 4- (9iMn> —$\ — — I 



d 2 d) , N l<90 [AiMa + a 2 M 2 + m 2 L 2 . 



Note that here we treat x~ coordinate special and replace all the d x - as iM, because this 
coordinate plays a distinguished role in the Schordinger holography. Later on we consider 
the parameter M to be either real or imaginary, from the observation that the coordinate 
x~ can be treated as an internal coordinate from the field theory point of view. 
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2.3 Two-point Correlation Function 



To find the solution of the equation (10), we use the Fourier decomposition as 

<P(u,x\y) = J e**T u {x + ) f^u) O M) , (11) 

where k is the momentum vector for the corresponding coordinates y. <po(u>, k) is introduced 
for the calculation of the correlation functions and is determined by the boundary condition 
with the normalization f^^ius) = 1- And T LU (x + ) is the kernel of integral transformation 
that convert u to x + , which is necessary for our time dependent setup. 



With this Fourier mode, the differential equation (10) decomposes into time dependent 
part and radial coordinate dependent part. The time dependent equation and its solution 
read 

^ + id+^T u =uT u — ► T u] (x + ) = c 1 exp- i ^ + (x + y J f L . (12) 

The radial dependent equation is given by 

u 2 f'^ + (2iMa - 3)uf'^ - (AiMa + a 2 M 2 + m 2 L 2 + u 2 k 2 ) = 2Muj u 2 j^ k , (13) 
where /' = d u f. An analytic solution is available in terms of Bessel functions as 

L % = u "' iaM (C2l»(qu) + c 3 K u (qu)) , (14) 



where I u and K u are Bessel functions with v = + L 2 m 2 and q = \/k 2 + 2Mu. 

For this solution to be well defined, we need to impose two different physical boundary 
conditions, one in deep inside the bulk and another at x + = 0, because we consider only 
x + > 0. To have a well defined field deep in the bulk, we choose K over /, 

L$ = c u 2 — M K u {qu) , (15) 
which is exponentially converging for large u, deep in the bulk. Thus the full solution is 

<P(u,x\y) = J |4^™V (^A 2 c K u (qu) . (16) 



iaM 



Thus we confirm that the aging scalar function is scale- invariant prefactor (^r-^J 2 times 
that of the Schrodinger background, advertised in [12] 



iaM 

4>Agmg{u,X + ,y) = ( —j- ) 0Schr6dinger(M,a; + ,y) . (17) 



We consider the x~ coordinate as the one for the internal manifold and restrict the system 
to be a sector with an eigenvalue M. The parameter M would be real if the x~ coordinate 
is periodic as in discrete light cone quantization. We assume that the x~ has boundaries, 
for example, at x~ = and take M to be a general complex number. This seems to accord 
with the the boundary at x + = 0. It would be natural to expect that the imaginary M 
would represent the dissipative behavior of the system. Let us consider the radial fall offs of 
the wave solution of the scalar field with the conformal dimension 2 — v. At the boundary 
u — y 0, the solution behaves as 

L,h ~ u 2 ~ u u~ iaM ~ u 2 ~ u+aM 'u- iaMR , with M M R ■ . (18) 

Note that we consider complex M = Mr + zMj for notational simplicity. Below, we will 
consider real M = Mr and imaginary M = iMj cases separately to make the distinction 
clear. 



It is worthwhile to pause and have some comments. First, the conformal dimension of 

itxM 

the scalar field </> is the same as in AdS in light-cone because the factor ( ) is scale 
invariant and thus does not contribute to the scaling dimension, which is observed in [12]. 
Second, as a result, the fall-off behavior of the scalar wave solution at the boundary changes 
for the case M = iMj. This explicitly demonstrates that the correct scaling dimension of 
some operators can not be read off just from the behavior of the radial wave solution in the 
time dependent holography. Third, this change of the radial wave solution can be viewed as 
"wave function renormalization" and should not be canceled by additional counter terms, as 
observed in [12j. Physical quantities would have appropriate time dependence to compensate 
this extra effect. This will be explicitly demonstrated by evaluating the correlation functions 
below. 

We further require the scalar wave solution to behave well as it approaches to x + — > 0, 
because x + = also serves as a boundary. This is similar to the boundary condition imposed 
deep inside the bulk along the w-direction, which leads physically reasonable results below. 



To check the behaviors, we calculate the wave solution, equation (16), near the boundary in 



coordinate space for the massless case m = with M = iMj, v = 2 and <po(u, k) — 1 as 



a.M 



Mjy 



772 



^(z+r^^exp j-^j . (19) 

This shows the typical behaviors of the wave solution in the coordinate space. One sees that 
the wave solution converges at x + = for Mj > due to the exponential factor and also 
converges at x + — > oo for + 1 > due to the polynomial factor. This demonstrates the 
convergence of the wave solution at the time boundaries for the imaginary M = iMj. From 
the fact that the asymptotic forms of the wave solution and two point correlation function 
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are restricted by aging symmetry, one expects to get similar exponential and polynomial 
factors for more general cases. 

We follow |34j [35J [36J [37J to compute the correlation functions by introducing a cutoff ub 
near the boundary and normalizing / — lj which fixes c = u~j 3 2+taM K~ x (qus) ■ Let us 

consider the on-shell action, which has two different contributions, the terms proportional 

to g uu and to g u ~ , as 

SIM = J d 3 x^r <P*(u,x + ,y) f^d u + iM^\ <P(u,x + ,y)\ UB . (20) 



This can be recast using the equation (16) as 

rlr.i' rlt.i iu(M*-M) 

dx + 9 (x+) ^L^e-^'-^ (ax + y^-^ 
2n 2w 

X J d2y J0^J ^^"^^(^^K^w.P.^Co;^!^ , (21) 



where 8(x + ) represents the existence of the physical boundary along the time direction as 
< x + < oo, and T is given by 

T(u,u',u>,fi,k) = ^fc )j? (u',&,u) (^d u + iM^j f u fi(u,k,u). (22) 

Note that the spatial integration along y can be done to give delta function 5 2 {k' — k). One 
can bring the u ±taM factors in / and /* together to cancel each other, which removes the 
second part in J 7 . Then using the relation 

d 

—u 2 K u (qu) = u{(2 - u)K v (qu) - quK^iqu)} , (23) 
au 

one can evaluate the u-dependent part at u = ub at the boundary by expanding in terms of 
small Ub to obtain the non-trivial contribution in J 7 as 

Note that the function J 7 is only function of u and k when it is evaluated at the boundary. 
Thus the time independent part is the same as previously considered cases in [1] [5][T2"]. 
But this is not the end of the story. To obtain the final form of the momentum correlation 
function, we need to evaluate the time dependent part. Here we consider two cases, M = Mr 
and M = iMi with real parameters Mr and Mi, separately. 



10 



2.3.1 M = M, 



R 



One can perform the x + integration in equation (21 ) to get the following results, because of 



the step function and vanishing exponent, M* — M = 0, 

,-■/ ; \ 5(oj' — uj) i , . 

Q(J -u) = ; r . 25 

Thus the momentum correlation function has two parts 

(0*(u f , k')0(u, k)) = -2{27c)- 3 5{k' - k) F{u B , u, k) Q(u' - u) , (26) 



where T is given in equation (24). This is the final form of the momentum space correlation 
functions. Note the extra term present due to our boundary condition with step function in 
time domain. 

Let us evaluate the coordinate correlation function aa3 

(<D*(x+, y 2 )<D(xt, ft)) 

<U <H '' ^ (1HZ ^.y 2 -ik-y, e -^.4^-4 2 (0*(u',k')0(u,k)) . (28) 



2tt (2tt) 2 2tt (2tt) 2 \xj 
For further calculations, we use the following integral 



*£ e -W-»)4g {u / - u) = l±*p?Il = e(xt) , (29) 



with the condition x^ > 0, and the inverse transform of q as 



1 -^-exp [iM R ^\ yx > , (30) 



7r2 1 -"i 1 +T(-i/) (x+ - xlY+» ^\ n 2(xZ-xj 



4 Note the time dependent factor (x + ) 1 2 in front of the boundary scalar operator is opposite to that 



of the corresponding scalar held in (16). The scaling dimension of the aging scalar operator is effectively 



different by iaM compared to that of the Schrodinger case 

[paging] = [^Schrodinger] + iaM , (27) 



as read off from the radial fall-offs in equation (17 1. This should be compensated by an appropriate time 



dependent factor. Another way to see this is the the relation 

d(po{l) 

where the time dependent factor would result in opposite way. This is the same for the case with M = iMj 
discussed below. 
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where we used {y 2 — Vi) 2 > and > 0. Then the coordinate space correlation 

function is 



da; d 2 /c _.-ir./-.-t._A-i „•,.,/■„+ /" da;' 



x+ / y 2vr (2tt) 2 y 2tt 

r(i-iz) l 3 m^ u e(x+)6(xt - x+) fxt\ 1 ^- f. M R (y2-yi) 2 \ , 01 , 



" r(i/)r(-i/) tt2^- 1 z 1 +^- 2i/ (x+ - x+) 2 +^ \xf J V 2(x+ - x+) 

This is our main result for the zero temperature correlation function with M = Mr. We 
would like to pause to comment: 

• We find that the time dependent polynomial part is the same as the result [12] when 
M = Mr. This is associated with the fact the the scaling dimension of a boundary 
scalar operator is opposite to that of the source of the corresponding scalar field. 

ictM q 

(0%x+,y 2 )0(x+,y 1 )) A ^ ng ^(^:j ' (0*(x+,y 2 )0(xf, £)> Schroder • (32) 

• It turns out that there exist only oscillating behaviors for M = Mr. Thus the same is 
true for real value for aM, because there exists only the combination aM. 

• For the massless case with u = 2, we have 

F(u B , u, k) = ~q A log(g 2 ) + ■ ■ ■ . (33) 
Thus the coordinate space correlation function is given by 
(0\xt,y2)0{xt,yx)) 

itfMl e(xt)e(xt-xt) fx+y^ f . M R (ih - i/i) 2 \ ^ 



/I V 2 1 / \ 1 / \ \ 2 1 

2.3.2 M = %Mi 

Let us evaluate the time integral first. We perform the Fourier transform for the time range 
< x + < 00, and define Q(w) = J dx + e iwx+ 6(x + ) ■ (ax + )~ aMl . Then 

Q{w) = a- aMl \w\- 1+aM 'T{l - aMi) f-zco S (^^)sign(u;) + s in(^^)^ . (35) 



where w = uj' — u. Note that these results are essentially the same as those evaluated by 
Laplace transform in, so called, s-domain. 
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Thus we get the momentum correlation function as 

(0*(u', k')0(u, k)) = -2{27r)- 3 5{k' - k) Q{J - u) F(u B , 00, k) . (36) 



Again, J 7 is given in equation (24), and q = \/k 2 + 2Mu. Here we can have k = due to the 
spatial translational invariance, but not for u because time translation symmetry is broken. 

Let us calculate the coordinate space correlation function as 

(0*(4,jft)0(4i)) (37) 
^S^** 8 '*"'*^ (a 2 xM) ^ (0*(u>', P)0{u>, k)) . 
After carrying out the integral for the u' first 

<^ e i(u/- U )zt Q^J _ u)= QU+\ a -a Ml ^- aMl ? (3g) 
Z7T 

the expression reduces to the previous case for the time independent part, whose calculation 
is done similarly. We get 

(cr(xt,jh)o(xt,jh)) 

= -20(4) (a 2 xfxt)~ a- aM '\xt\~ aMl / ^^e^^e^-^ F{u B , u, k) 

J 2rc (zti) 1 

otMj 

T{l-v) L 3 M} +U e(x+)e(x+-x+) fx+Y^ ( M^-y^ 

' eX P 1 r>/ + +\ I ■ 



T(u)T(-u) 7 r2^ 1 u B - 2u (x+ - xf) 2+u \xf J \ 2(x 2 - x{ , , 

This is the result we are looking for. This confirms that the correlation function has dissi- 
pative behavior for M = iMj. Thus 

• We find that the time dependent polynomial part is identical to the result [12] when 
M = Mr. This is associated with the fact the the scaling dimension of a boundary 
scalar operator is opposite to that of the source of the corresponding scalar field. 

aMj 

(0*(x+,y 2 )0{xf,y 1 )} A ^ ng = (^j ' (0*(x+ y 2 )0{x+, yi)>Schrodi»g e r ■ (40) 

• There exist only the combination aM, and thus the result is true for imaginary aM. 

• For the massless case with v = 2, we have the coordinate space correlation function 
similar to the M = Mr 

olMj 

{O (xj,y2)0{xj,y 1 )} = - + + — exp — -— ^ 

7T [Xj — XJY +V \xj J \ 2(^2 — X{ , 



(41) 
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Figure 1: Plot of the equation (42) with v 
and vertical axes are Ax + 



0, aMj = 1, aMft = 0. The horizontal 
and two-time correlation function C = C(x~2,xt), 



respectively. The four plots, black, red, green and blue (from left to right) are for the 
particular waiting time x^ = 25, 50, 100, and 200, respectively. They show that the older 
is the system the slower it relaxes. This plot is to be compared with that of the numerical 
results in [TO"] . 



2.4 Two-time Correlators 

One of the important physical observables considered in the literature (see e.g. |9][3E]) is 
two-time correlation functions at the same space positions, which we refer to C(a;^,x^) 

C(xt,xt) = (0*(xt,y 2 )0(xt,m))y 1= y 2 = (xf)- x » f c (^) , (42) 



X 



where fc (x) = x~ x °/ z , and z is the dynamical exponent, xf and are called a waiting 
time and a response time, respectively. This scaling behavior is expected to apply for the 
"aging regime": xf,x^ 3> x^i and x\ — x\ ^> x^, where x^ is a microscopic time scale of 
a given system. 

Our dual field theory possesses aging invariance, inherited from the bulk gravity and hence 
exhibits two important characteristic features of aging system, namely the time dependence 
and the existence of dynamical scaling. Furthermore, for Mr = and for x\ 3> xf , which is 
physically interesting aging regime, our two-time correlation function decays slowly following 



power law, like (42), with 



, aMi s , 

A 6 = v + 2, A c = z \v + 2 + — - , (43) 
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where z = 2 for our case. These two parameters A& and Ac are our physical exponents asso- 
ciated with the two-time correlation function. The dissipating behaviors of the correlation 
function are plotted in figure [T] with several different waiting time. We clearly see the typical 
aging behavior: older systems relax in a slower manner than younger systems, and thus the 
correlation between two different time x^ and x\ has more correlation for the bigger waiting 
time x\ in our correlation function C{x\, xf). 



2.5 Comments on Schrodinger Background 

In this section, we would like to briefly comment on the aging construction on the Schrodinger 
background [12] with our approach, especially with the complex M and without complexi- 
fying the spacetime geometry. 

Schrodinger background can be embedded in string theory in a nice way using null Melvin 
twist [19] [20]. Previous sections already confirmed that the results of the scalar wave solu- 
tion and the two-point correlation function are the same for both cases even though these 
two approaches are rather different. This is mainly because the Schrodinger symmetry is 
large enough to restrict the physical properties of two point correlation functions. It is 
still worthwhile to revisit the aging construction for the Schrodinger background with our 
approach. 

Let us consider the following metric with the change x + — > t, u — > z and x~ — > £ 

ds 2 = ^ (dz 2 - \dt 2 - Idtdi + dy 2 \ + ^ f— dzdt - - dtA , (44) 

Z y Z J Z V Z t J 

where extra matter contents, which are not explicitly written here, are necessary to support 
this background. Following the previous section, we couple a scalar field to this background. 
We find the identical radial differential equation as the one in AdS in light-cone if we change, 
m 2 m' 2 = m 2 + M 2 /L 2 [5]. 

For the complex parameter M = Mn + iMj, the form of the corresponding solution is the 



same as dl6h with q = VW + 2Mu and v' = a/4 + L 2 m 2 + M 2 . These parameters were used 
in the Schrodinger background [5] [1] . Again we confirm that the aging scalar wave solution 

iaM 

is given by the scale-invariant factor, (ax + /u ) 2 , times the Schrodinger wave solution as 



advertised in [12]. This is given in (17). With this in hand, it is straightforward to evaluate 



the two-point correlation function, which is also the same as that of the AdS in light-cone and 



is given in equations (31) for M = Mr and (39) for M = iMj. Thus, at zero temperature 



the physical properties of Aging in light-cone are the same as those of the Aging background 
for the operators with the same scaling dimensions, v = v' . The scaling dimension of the 
operator is complex in general for complex M. If we consider either real M = Mr or pure 
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imaginary M = iMj, which are physically interesting cases, the scaling dimension is still 
real. It is interesting to notice that the scaling dimension — 2 — v' decreases for the real 
M = Mr, while it increases for imaginary M = iMj. 

All the other discussions and physical properties remain unchanged compared to the 
Aging in light-cone at zero temperature. Interestingly enough, as we will see later, these do 
not persist at the finite temperature. 



2.6 Comments on Dual Field Theory 

In the context of the holographic condensed matter application, the field theory dual of the 
gravity description is not explicitly known in most cases. For the case of the Schrodinger 
geometry, the situation is much better because the construction are explicitly known as 
Discrete Light-Cone Quantization (DLCQ) for AdS in light-cone and null Melvin twist for 
the Schrodinger background. The basic field theory dual of AdS in light-cone is DLCQ of 
the A/" = 4 Super Yang-Mills theory proposed in [39J [40J as explained in [TS] . It is further 
required to project the conserved particle number M to a single sector, which is related to 
the isometry of the gravity backgrounds and has distinguished role [31 El El [26]. We are 
required to focus on a single sector of the particle number. For further developments along 
the Schrodinger background with null Melvin twist, see \\J)\. In this section we would like 
to comment on some properties directly related to the time dependent generalizations of the 
setup, which apply for both AdS in light-cone and Schrodinger backgrounds. 



The differential equation (10), we consider at the end of section 2.2 has an explicit time 
dependence. Thus the dual boundary field theory is expected to have not Schrodinger but 
aging symmetries. From the observation that the equation ( JlOp factorize into u— dependent 
and u — independent parts, we factorize = ip(x + ,y) f(u) and find that the boundary wave 
function tp satisfies the following differential equation 

2M (i^- + cp + V 2 cp - v 2 cp = , (45) 



dx + 2x + r 

where v 2 is the eigenvalue of the radial differential equation. The term v 2 <p can be interpreted 
as the one coming from the variation of the Landau-Ginzburg potential, i.e. ^ = v 2 cp. 
One may generate the general Landau-Ginzburg potential by considering the self-interacting 
scalar fields in bulk. It is clear that our geometric realization of aging symmetry actually 
constrains the nature of the time dependence of the equation of the dual field theory as the 
form 

aM 2 . , 

— V . (46) 
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It is interesting to note that if we can add noise contributions the differential equation 
becomes the Langevin equation 



It is not difficult to solve the differential equation (45) and read off the time dependent 



part of the solution. Compared to the a = case we get 

I ia.M 

ip — > (x + ) 2 ip . (47) 



Thus with aging symmetry, the boundary scalar field will acquire the explicit time-dependence 
of this form. Note that we can also infer the same conclusion from the local transformation 
(J3j) which contains essential effects to the bulk side of the story. Furthermore, the parameter 
v contains other effects due to the renormalization of the radial coordinate. It turns out that, 
for real M, there is no net effects of this radial wave function renormalization on physical 



observables such as the two-point correlation function as we see in (17) and (31)(39). 



3 Finite Temperature Aging Holography 



We start with the Einstein- Hilbert action with cosmological constant, which admits the 
planar black hole solution 



+ 



^ \ 2 
L 
L 



1-h 
4b 2 



-dr 2 



(dx 



[1 + h)dx + dx~ + (1 - h)b 2 (dx- f + 



h=l- 



H 



(48) 



We would like to exploit the idea that the aging metric is locally Schrodinger, and thus 
we can get the Aging black hole from the Schrodinger black hole using the coordinate trans- 
formation as 



x — y x + —In (r 2 x + ) 



(49) 



The black hole metric in the simplest form can be written as 



ds< 



L 



-^-b~ 2 dx +2 + (1 - h)b 2 ( dx~ + -dr- , 
1-h v ; V r V 2x+ 2b 2 1 - h 



L 



+(-| h~ l dr 2 
r 



+ \ + ^ ) dx + ) +dy 2 



(50) 



We would like to consider this finite temperature metric to construct the correlation func- 
tions. 
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3.1 Equation of Motion 



We consider a probe scalar on the Aging black hole background and try to evaluate the 
two-point correlator similar to the previous section. The equation of motion is given by 

' -A W^dsTdut) - m 2 (j) = . (51) 



V-9 

The metric depends only on (x + ,r) and the equations of motion can be expressed as 

^(4 + h)d r <f> - -cL0 - ^^y 2 — d-4> + - "V = • (52) 

As was done in the previous section, we replace the derivative of the x~ coordinate as 
d x - = iM. The mode expansion becomes of the form, 

<f>(r,x + ,y) = J ^^-^exp-^0 fc (r,x + )0 o (w,fc) , (53) 

where (p k (r,x + ) = T u (x + )f u ^(r). 

The equation of motion in r coordinate is complicated. It is slightly better to go to a 

r 2 

different coordinate using z — -f-. Then the differential equation is 
l + h M DT ul (l-h)b 2 D 2 T u 



4zhr 2 H V h "' T„ h T L 

_ gj (-iaMh -h + 2) fls ( + a 2 M 2 h + m 2 L 2 fc 2 - 

4 g ^ / WfS ^ 4^ + 4^ ; • 1 j 

where /' = c^/, and we used 

d olM o . . 

One qualitatively different feature from zero temperature case is that this finite temperature 
version of the differential equation does have a second order time derivative term. Yet as we 
see in the first line of equation (54), this does not generate different dynamics along time 
direction compared to the zero temperature case because the time dependent parts of the 
differential equation, equation (54) factorize in a special way. 

The time dependent part can be solved as follows: 

'aM 



2x- 



+ id + ) T w = uTu — ► T w (x + ) = Cl exp-^ + (x+)^ . (56) 



Thus, we would like to stress that, the time dependent part of the scalar wave solution at 
finite temperature is exactly the same as that of the zero temperature case. Anticipating the 
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appearance of the scaling invariant combination, x + /z, we factorize as f u %(z) = z ^ g u %{z) 
and obtain the radial equation: 

„ 1 + z 2 , m 2 L 2 „ w 2 z „ q 2 
9 m ^9 ~ t^tm ^9 + 7: ^9 m = , 



2\ 



z{\-z 2 Y Az 2 (l-z 2 Y (X-z 2 ) 2 * z(l-z 

M/(2b 2 )-u 2 2Mu + k 2 

where w = — — — ^ , a 2 = — — Im .„ , (57) 

2ttT ' H {2ixbT) 2 ' v ; 

where used r# = nL 2 bT. The equation includes the term which has the Schrodinger invari- 
ant combination 2Mu + k 2 , explicitly, and the term which is partially broken due to the 
other combination w ~ Mj (2b) — boo at finite temperature. The term proportional to w is 
dominant at the horizon z — 1 and important for the calculation of the correlation functions. 
Interestingly enough, the combination of chemical potential ^ and conserved charge M of 
x~ direction of the time-independent Schrodinger black hole shows up with the energy of the 
probe scalar oj. In what follows, we consider the low energy regime in which the magnitude 
of the combination ^ is bigger than u, which corresponds to w > 0. 

Note that this equation is nothing but the radial equation for the AdS planar black hole in 
light-cone. Thus we explicitly check that the radial part of the time independent aging wave 
solution is the same as that of the AdS in light-cone case. This has an important implication 
that the time independent part of the aging wave solution and correlation functions are those 
of the AdS in light-cone at least in momentum space with some modification, which can be 
identified and thus isolated. 



The full solution is 



iaM 

d k duj ■ - ' ^ 2 



<K*,* + ,if) = J (|^i^ e *^ + (^rj c v^) , (58) 

where the solution has the desired scaling form with an appropriate normalization for the 
correlator calculation. We concentrate on the momentum space correlation functions from 
now on. 



3.2 Finite Temperature Correlation Functions 

In this section, we would like to evaluate the scalar correlation function of the Schrodinger 
geometry via AdS in light-cone with finite parameters M, oo and k, which is fairly non- 
trivial even in this simple setup. First, we consider the time dependence of the coordinate 
space correlation functions and outline the general properties of the correlation functions at 
finite temperature. After that we concentrate on the momentum space correlation functions 
starting with analytic approaches in general. And then we evaluate the correlation functions 
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analytically for some special cases and obtain the expression for the shear viscosity, which is 
presented in appendix [Cj Finally we get the full features with numerical studies. From the 
numerical study, we observe two distinct peaks in momentum space correlation functions, 
one broad peak at the (ou, k) = (0, 0) and another at finite values of to and k, which are 
explained in detail at the end of this section. 



3.2.1 General Features 



In section (3.1), we have some important observations regarding the time dependence of the 



scalar equation of motion. Combining with the experience from the zero temperature case 



in section (2.3), we would like to see the general time dependence of the coordinate space 



correlation functions. 



Following the steps done in section (2.3), we rewrite the solution (58) as 

d 2 k duj 



where / Az) = cz~ % ~ s ^ L g u) t{z). Then the on-shell action has the form 



(59) 



(60) 



This can be recast using the equation (59) as 
, s duj' dcj 



dx 



[x 



-i(to'— L})x+ 



x / d A y 



2tt 2tt 
d 2 k' 



in(M -M) 
2 



d 2 k 



(2tt) 2 J (2tt) 2 

where 8(x + ) represents the physical region as < x + < oo, and J 7 is given by 
T{z, to', 9, k) = yf^f*. 9, z) (g zz d z + iM g z ~) f w% (u, k, z). 



(61) 



(62) 

Note that the spatial integration along y can be done to give delta function 5 2 (k' — k). One 
can bring the z ±l9 ^ factors in / and /* together to cancel each other, which removes the 
second part in J 7 . Then time independent radial part J 7 of momentum correlation function 
is given by 



J-(z, cu', to, k', k) 



6 4 7T 4 L 5 T 4 

2^3 



Az 2 h 
~1? 



-d z 



(63) 



It turns out that evaluating this correlation function is highly nontrivial because the differ- 
ential equation is not analytically tractable for the AdS in light-cone except some special 
limits, while it is even more complicated for the Schrodinger background. 
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Before explaining the time dependent part of the correlation function, we would like to 



comment on our choice of the incoming boundary condition for the equation (57), which is 
not clear due to the contribution of M in the w. We choose our incoming boundary condition 

as 



[\-z)™/ 2 <x e -M* + +**) (64) 



where z* = 4 T ■ We come back to this radial dependent part of the correlation function 
with several different approaches below. 

As discussed in the section [2j this is not the end of the story. The time dependent parts 
turn out to be the same as those of the zero temperature case because their time dependences 



in the wave solution (58) are the same. This is explained below equation (55). Thus the 



time dependent parts of the correlation function can be explicitly computed and are given 



by Q in equations (25) and (|35j). Explicitly, we get 

(0*(u', k')0(u, k)) = -2{27r)- 3 5{k' - k) F{u B , u, k) Q{J - u) , (65) 



where T is given in equations (63). By performing the u' integration, the coordinate space 
correlation function is give by 

(0*(xl,y 2 )0(xf,y 1 ))=9(xZ) (^J " (0*(x+, y 2 )0{xt, Schrodinger - (66) 

Where M can be either M = Mr or M = iMj. The latter case is related to the dissipative 
case. 

It is not straight forward to evaluate the general momentum space correlation functions 
analytically, not to mention the coordinate space correlation functions. Thus we would like 
to concentrate on the momentum space correlation function for the rest of the section. 

3.2.2 Analytic Approaches 

Let us outline the analytic approach we considered and summarize our observations. 



The differential equation (57) has 4 regular singular points, which, in general, can not be 



solved with systematic approach. If we take 



g(z) = z 1 ¥{z-l)-% 1 {z + l)- J FG(z) . (67) 



The differential equation (57) can be cast into the Heun's equation 



G "« + (? + + ITl) G ' {z) + sts-iHsV *" = • (68) 

where 4af3 = v 2 - 2 + 2(7^ + 8e + ey) , -2q = q 2 + 7<5 + 7e , (69) 
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with the special values of 

7 = 1 ± v , 8 = 1 ± iw , e = 1 ± w . (70) 

The general solutions of this equation are known as the Heun functions Hl(—1, q; a, {3, 7, 5; z) . 
In general, we have 8 different combination of the parameters 7, 8, out of which only 4 
solutions are physical due to the incoming boundary condition we would like to impose at 
z = 1. That corresponds to the choice 5 = 1 + iw due to our choice of w parameter given in 



equation (57). For special values of parameters, the Heun functions might provide analytic 
solutions. 



We observe that our differential equation (57) is similar to the equation (3.4) of [42J, 
even though they have very different setup. We can rewrite the differential equation (57) in 
a slightly different way to make the comparison straightforward 



~ 9 " - w^)~ g ' - U 2 (i - ^ + 41 - ^ )~ g=0 - (71) 

where Y 2 = P+{ ^ ] T t hui) \ Q 2 = w^f • Note that, in general, Y 2 > Q 2 , where the equality 
holds when M = 2ub 2 . One of the crucial difference between this equation and (3.4) of |42] 
lies in the last term. In our case, Y 2 has always larger contribution than Q 2 , and the term 
proportional to Y 2 is multiplies by z 2 . On the other hand, the corresponding term in (3.4) of 
[42] can be neglected to give a constant numerator, which make possible to proceed further. 
Due to the difference, we can not use the approach adapted in [42] in a reliable way except 
for a very special case, which we don't pursue further here. 

In subsequent sections, we consider high temperature limit, called hydrodynamic limit, 
and low temperature limit in some details. And then we evaluate the differential equation 
using numerical method to see the full features of the correlation functions. 

3.2.3 High Temperature Correlation Functions 

For the high temperature limit, we have small parameters w and q 2 , which can be used 



as expansion parameters in equation (57). Following [34], we can derive similar results 
[14] [T6] [26] . For m = case, we calculate T using incoming boundary condition at the 
horizon 



T = 26VL 3 T 4 + (iw + q 2 ) + O (w 2 , wq 2 , q 4 ) J 



(72) 



where we used 



, ,(.) = (1 - zr 12 ( 1 - (iZ + q 2 ) In ^ + ■ ■ ■ ) • (73) 



2 V 2 
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Note that the form of the first non-trivial contributions are the same as the results of [34] 
even though the differential equations are different. Yet, the physical implication seems to 
be more complicated. In general we can not claim to have q = for the zero momentum 
case k = because of the parameter M. The first term diverges for the strict ub — > limit, 
while the second term gives finite contributions as 

Where we put the result with M = Mr + %Mj for notational simplicity. 

Let us concentrate on the case q = 0. This explicitly means that 2M^u + k 2 = and 
Mi = 0. Thus the momentum correlator has an extra contribution in imaginary part of the 
correlator compared to relativistic case [31] 

Re T = , Im 7 = 26 4 vr 4 L 3 T 4 Mr ~ ^ . (76) 

Anb 2 T 

If one considers the case of zero spatial momenta k = and thus M R = 0, we get 

Im J= = -b\ 3 L 3 T 3 u . (77) 
For the application to shear viscosity, see appendix [Cj 

3.2.4 Low Temperature Correlation Functions 



With a change of the function g(z) = J z — - g{z), we get the differential equation for g(z) 



g"={w 2 F + H)g, (78) 
p __ (1 — z 2 )s 2 - z 2 H= -3 + 6z 2 + z 4 -m 2 L 2 (l-z 2 ) 
z(l-z 2 ) 2 ' ' " Az 2 (l-z 2 ) 2 ' 

where, in the low temperature regime, q 2 ,w 2 ^> 1 and s 2 = ^ = (M^/^^-buj) 2 ~ ®(^)- 
Then the term proportional to w 2 dominates the potential, and the solution can be found by 
the WKB approximation following [31]. We concentrate on the case Mj = in this section. 



The Schrodinger equation (78) has a potential, which is positive when < z < zq and 



negative when z < z < 1, where 



= 2M ^ + fca (80) 

v/ 1 + s 2 V (M fl /(26) + 6w) 2 + /c 2 
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Thus the solution of the equation (78) decays exponentially in the interval < z < zq and 
oscillates in the interval Zq < z < 1. Physically the particle has to tunnel from z = to 
z = Zq before it can reach the horizon z — 1. The imaginary part of G R is proportional to 
the tunneling probability 



Im G H ~ exp ^-2|w| J dz^F{z)J , (81) 
and this integral is easy to evaluate for the following two limits : 



For small sj^jwe get f*° dz^F\z) = ~ r ^\y^ (s 2 ) 3/4 , and 

t rR ( 0.556 (2M R co + k 2 )^ \ 

ImG = eX »{-^- (M R/m - b ^ ) ' (83) 

For large s, we get f*° dzy/F(z) = r /3 \ an d thus 



,2M R u + k 2 

Im G R = exp ( -0.835 1 — | . (84) 



Thus, at low temperature, we obtain the generic behavior of the imaginary part of the 
correlation functions as an exponentially decaying one, e~ A ^ T . Specifically, the result for 
large q, where Im G R ~ exp (— q), is further confirmed by our numerical study in below. 



3.3 Numerical Results 

In this section we would like to present the numerical results of the correlation functions by 



solving the differential equation (57). For massless m 2 = case, the series solutions of the 
boundary and horizon are presented in appendix [Bj Massive case has also similar solutions 
with an additional parameter m 2 L 2 , whose behaviors are briefly explained below. 



5 One might wonder that this condition invalidates the assumption q 2 ,w 2 3> 1. For w > q, one can 
evaluate the following differential equation 

9 "=U*F + H) 9 , F = (1 - f } - J 2 , H = H. (82) 
v ' z (1 — z z ) 

where s 2 = 1/s 2 > 1. Then the term proportional to q 2 dominates the potential, and the solution can be 
found in a similar manner. 
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to 



Figure 2: Real and Imaginary parts of the scalar correlation function in terms of the pa- 
rameters k with uj = for left plot and u with k = for right plot for the fixed values 
M = 10, b = 1 and T = 1. There exist well defined peaks in imaginary part of the correla- 
tion function for both cases. 



3.3.1 Description in terms of uj and A 

For fixed M — 10, b = 1 and T = 1, we present plots for the real and imaginary parts of the 
two-point correlation function in figure [2] The left plot is for the fixed uj — 0. There exist 
a sharp peak at finite momentum and also a much broader peak at zero momentum k = 
for zero frequency. As we increase u to the positive value, the peak moves toward to the 
k = and the peak becomes bigger and narrower, while it becomes smaller and broader as 
we decrease u to the negative value. The right plot is for the fixed k = 0. There also exists 
a sharp peak at finite value of u for zero momentum. As k increase, the sharp peak moves 
toward negative u and becomes smaller and broader. Thus, overall, the sharp peak becomes 
broader and smaller as we increase k and decrease u. 

These behaviors can be checked in the figure [3j Note that we exclude the region which 
gives imaginary q. There is another much broader peak at (u>, k) = (0,0), which can be 
identified as the point q = 0. This peak becomes smaller and broader along the line identified 
by q ~ 2Mw + k 2 = 0. This peak can be clearly seen in an alternative description below, 



section 3.3.2 The above mentioned sharper peak seems to follow the broader peak line with 



some fixed distance, which can be described as 2Mui + k 2 = This peak is identified as 
the case where w = 0. Thus we confirm that these two peaks are related to the case when 



each of the last two terms in (57) has vanishing contributions. 

Let us comment about novel features of the scalar correlation functions for the Schrodinger 
geometry. First, we clearly have two peaks, of which the higher peak is distinctively clear 
and directly related to the unique property of the Schrodinger geometry, the isometry M and 
associated chemical potential ^ of the spectator direction. This, apparently, is analogous 
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Figure 3: Imaginary part of the scalar correlation function in terms of the parameters u, k 
for the fixed values M — 10, b — 1 and T — 1. They are even functions of k and we take the 
parameter range as —10 < u < 5, < k < 15. The plots are only drawn for the parameter 
range 2Mu + k 2 > 0. 



to the quasi-particle peak at the Fermi surface in the case of spinor. Second, the shape of 
our peaks in the (u, k) space follows a curve 2Mu> + k 2 = const, due to the typical energy 
dispersion relation in non-relativistic Schrodinger type, rather than relativistic ones which 
would be linear in the same parameters. 



3.3.2 Description in terms of q and w 

While it is more physical to parametrize the physical properties of the correlation function 
in terms of u, k, some of its properties are more transparent in terms of the parameters, w 
and q. Here we add some comments on them. 

The typical shape of the correlator is depicted in figure |4j We would like to identify the 
corresponding peaks we observed in section 3.3.1| One of the peaks is located in q = 0, 



which is matched to the broad peak located in u = k = 0. This is not exactly located at 
the origin but located in slightly non-zero value of the parameter w = ^ 2T , due to the 
non-trivial parameter M for the Schrodinger holography. The other peak is very sharp in 
terms of the parameter q and spread out in w. It is located at finite q 2 = M ^ b ^y. , we can 
identify this point near w = 0. As we increase the value of w, the peak gets smaller and 
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broader, but the peak location does not change much in q. 

Let us comment about the behavior of correlation function as we change the bare mass 
of the scalar m 2 L 2 , which correspond to the change of the conformal dimension of the 
corresponding Field theory operator. For the range —4 < m 2 L 2 < —3, there is a single 
peak at q = and at small w, which spread out as w increases. As we increase m 2 L 2 , the 
single peak moves to some positive value of q without being separated from the peak at — q. 
Further increase of the bare mass brings the peak to the origin around at m 2 L 2 = 0. For 
positive mass range, m 2 L 2 > 1, there exist two separate peaks which are similar to the zero 
mass case in figure |4j while the peaks are small compared to the zero mass case. With the 
numerical analysis, we find that massless m 2 = case is rather special and has a well defined 
isolated peaks at q = and w = and spreading as we change the parameters. 

Finally we would like to comment on the correlation function at large q and large w 
with q > w. The figure [5] represents a typical behavior in the region. We can check the 
exponentially decaying behavior we obtained with analytic analysis in equation (84). 
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Figure 5: A typical log plot of the imaginary part for w = 3. 



3.4 Comments on Schrodinger Background 



In this section we consider the aging holography from the generalized Schrodinger background 
[14|[15j[16j generated through the null Melvin twist. The action is given by 

I /' ., . /,.. 4 
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16ttG 

where A and F are a massive vector field and the associated field strength, and is a dilaton. 
The black hole geometry associated with the action was found as 

1 



ds' 
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where 



h = 1 



H 



k = l + 



b 2 r% 



The massive vector field and dilaton also have the non-trivial configurations as 
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Following the previous section, we use the local coordinate transformation (49) to gen- 
erate Aging black hole inherited from the Schrodinger background as 



ds z 
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with the massive vector field 



.4 



r 
~k: 



2 r 



1 + h 



ab 2 rjj 
2x+r* 



dx + ¥-dx~ -^-dr 



We can check that the metric reduces to that of [12] when — > and also to that of 
|14j|15j[T6] when a — > by construction. One can make sure that these configurations 
satisfy the Einstein equation. 



We consider a probe free scalar on this background following the section 3.1 with the 



same equation of motion (51). Using the mode expansion 

duj d 2 k 



2tt (2 



7T 



exp 



S9) 



we can separate the x + dependent part of the differential equation with the same operator 



D given in equation (55). The solution of the time dependent part is exactly the same as 



before, T UJ (x + ) = ciexp lulx (x + ) ™ . This is one of the main observation for the aging 
backgrounds compared to the Aging in light-cone. 



Upon changing z — -| and / 
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(90) 



with the same temperature relation rs = TtbT and with L = 1 for the notational convenience. 
Again, this differential equation is identical to that of the Schrodinger backgrounds, confirm- 



ing the general relations (|17|), (|32|) and (|40|) in the Schrodinger black hole backgrounds. From 
these observations, we 

time dependent factor 



these observations, we argue that the local coordinate transformation (|3|) brings the same 

i a. Al p 



for M = Mr given in equation (32) and 



given 



in equation (40) to the correlation functions for the holographic aging phenomena. 



This equation (90) is the same as the corresponding equation (57) of the AdS in light-cone 



except the following change 



m 



m 2 k(z) 1/3 + M 2 , 



(91) 



which is another main observation of this section. We would like to comment about the 
consequences of this change without detailed investigation of this differential equation. 

• In general, k depends on the radial coordinate, which is crucially different from the 



case of AdS in light-cone. In there, the corresponding term in equation (57) is constant 
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and thus it is possible to consider the case m 2 = consistently. On the other hand, 
it is not possible to ignore this term for the aging Schrodinger case in general with 
non-zero M. This seems to be the first non-trivial evidence that there exist differences 
between the two different geometric realizations of Schrodinger holography. 

k is non-singular and monotonically increasing function from k = 1 at the boundary 
to k — 1 + ii 2 b 4 T 2 at the horizon. 

This is the finite temperature generalization of the difference between the AdS in light- 
cone and the Schrodinger backgrounds, m 2 — > m 2 + M 2 , which remains to be true at the 
boundary even at finite temperature. Thus we conclude that the differential equation 
of the Aging background is the same as that of the AdS in light-cone at the asymptotic 
region. 



It seems that the differential equation (90) is similar to the massive case of Aging in light- 



cone given in (57). It remains to be seen the differences in detail with further investigations, 



which we postpone for the future. 



4 Conclusion and Outlook 

In this paper we studied the time dependent aging systems by considering their dual gravity 
solutions with relevant symmetry, namely aging symmetry. We mainly focused on AdS in 
light-cone for the field theory with Schrodinger symmetry and its aging generalizations to 
incorporate the aging phenomena. Using these geometries, we study the two-point correlation 
functions of the dual scalar fields in aging system at zero temperature as well as those in 
Schrodinger system at finite temperature. At zero temperature, the two-time correlation 
function of the dual scalar fields exhibit all the relevant features of the typical aging system. 
Among others, it shows that the older is the system the slower it relaxes. At finite, but low, 
temperature with Schrodinger isometry, the analysis of momentum space correlation function 
exhibits the exponentially decaying behavior as a function of the inverse temperature. 

One of the important features in the Schrodinger holography and also in aging holography 
is the existence of the spectator direction in the bulk geometry. One usually chooses an 
eigenvalue M in that x~ direction, which appears as an extra free parameter in the boundary 
conformal field theory. In order to deal with the time dependent, open and dissipative system, 
it is tempting to relax the nature of this spectator or, more appropriately, internal direction 
and allow M to be a complex eigenvalue. Indeed, we could see that all the relevant features of 
the two-time correlation function in boundary aging system are revealed from the imaginary 
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part of M. 

Another important point in this paper is the clear differences we found between Aging in 
light-cone and Aging backgrounds at finite temperature. They are the conformal dimensions 
of the scalar operators and the mass dependence of the differential equations, which can be 
understood simply from the differences in the effective mass, m 2 in the former and m' 2 = 
m 2 fc(r) 1 / 3 + M 2 for the latter. For AdS in light-cone, we evaluate two-point correlation 
functions in momentum space for massless scalar case, m — 0, with finite M using numerical 
method and analytic approaches for two special cases, low and high temperature limits. It 
seems to be possible to connect this case to read off the viscosity in AdS in light-cone with 
finite M for general setup, while that was not possible in Schrodinger background due to the 
presence of w! in the fluctuation of the gravity modes. Further systematic investigations are 
required to see the details, while preliminary analysis is presented in [Cj Thus AdS in light- 
cone and its generalization to aging seems to be more attractive in their physical applications 
at finite temperature on top of the simplicity and well-defined boundary structure. 

Surely, we still need to clarify many issues to use the AdS/CFT correspondence for the 
time dependent field theory. Further extensive study of the aging system in this paper may 
shed some lights on the study in this direction. One particular example is the thermodynamic 
properties of time dependent backgrounds, which are readily calculated in this setup and 
appropriate applications can be clarified in a controllable setting due to the large isometry 
and simple setup [35] . 
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A High Temperature Correlator : Next Order 



In section 3.2.3 we evaluate the scalar two-point correlation functions for the first non-trivial 
order. We evaluate the next order with the same setup. Here we summarize only the results. 
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(93) 



Where PL is PolyLog. Following the standard prescription, we normalize the wave solution 
to be unity at z = and use the incoming boundary condition at z — 1. 

From the first two non-trivial orders of the expansion, we can calculate the momentum 
correlator as 



jr Q 2 

— - iw + q 2 + q 4 + (2iwq 2 + w 2 ) ln[2] - q 4 ln[4z B ] + 
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z B 



(94) 



Note that this is similar to the relativistic results jM] if we identify the parameters our w 
and q as relativistic frequency and momentum, respectively. But the physical meaning of 
these two parameters are very different and we develop several imaginary contributions by 



introducing complex M as we explained in section |3.2.3| In particular, there is another 
non-trivial contribution w 2 In [2] in the next order in addition to — iw even in the case q = 0. 
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B Series Solutions with m = 



Here we solve the differential equation (57) at the boundary and at the horizon with appro- 



priate expansions in each regions. These series solutions are used in the section [3T3] It is also 
interesting to consider general cases with non-zero mass m 2 , we would like to concentrate 
on m = in the differential equation (57). These series solutions are used to evaluate the 
radial part J 7 of the correlator as 
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At the boundary, we found a series solution upto 6th order of the expansion as 
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We explicitly list here only the first four terms as 
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This solution has two independent parameters, A and B, which can be fixed by the solving the 
differential equation at the horizon using incoming boundary condition and by interpolating 
the solution to the boundary. 

At the horizon we have a similar series solution upto 6th order of the expansion 
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with incoming boundary condition at the horizon as explained in the main body. The first 
four coefficients are 
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These two solutions, near boundary and near horizon solutions, can be connected numerically. 
And the two-point correlation function is given by 

J= B . . 

(100) 



46% 4 L 3 T 4 A 
in terms of various parameters we are interested in. 



C Shear Viscosity 

The finite temperature hydrodynamic correlation functions for the off-diagonal component 
of the energy-momentum tensor (T y f 2 T y ^ 2 ) were evaluated for the Schrodinger background 
in [H][T6], while that of the AdS in light-cone was done in [26J. It is interesting to recall 
that the viscosity-entropy ratio for the strongly coupled non-relativistic plasma also satisfies 
the universal bound, ^ = ^ [H][l6j[26j for u/T < 1 with vanishing parameters including 
M as confirmed in section 



3.2.3 



If we demand M — k — 0, the calculations are simple and 
effectively identical to those with the wave equation for a massless scalar field of AdSs black 
hole. 

Let us revisit the earlier treatments to evaluate the shear viscosity concentrating on 
Schrodinger geometry, which essentially used the scalar correlation functions, as advertised in 
|44j . While it is straightforward to identify the correlation functions of the energy- momentum 
tensor to those of scalar operators for the relativistic case with the vanishing spatial momenta, 
there exist subtleties due to the internal direction in the non-relativistic Schrodinger case. 
It is argued in [TJ] that M = is required to evaluate the shear viscosity in the Schrodinger 
background. The scalar operator has the conformal dimension A = 2 + \/A + L 2 m 2 + M 2 in 
the Schrodinger backgrounds, while the gravity modes have the conformal dimension A = 4. 
Thus it is required to put M = m = to calculate the shear viscosity using the scalar 
operator, and thus restricted to the sector, M = 0. 

The story is different for the AdS in light-cone. As we mentioned in several places, we 
concentrate on the fixed sector of M and we would like to consider the fluctuation equation 
of the off-diagonal component as in [Tj 



which satisfy the following differential equation 

+ 4 + ft s , + (M' + 4frV)(l-ft)-4^AMl + % = „ , m 
rh 4b 2 r 4 h 2 

Thus we are supposed to solve this differential equation directly to get the correlation func- 
tions related to the metric fluctuations. This equation is nothing but the scalar differential 



34 



equation given in (57) with m = 0, expressed in the ^-coordinate. Now a crucial differ- 
ence comes into play. Note that the conformal dimension of a scalar operator is given by 
A = 2 + v / 4+IV. Thus we can just restrict ourselves to the case m = 0, which is already 
taken care of in (57). This made it possible to read off the viscosity, for general sector M, 
readily from the boundary two-point correlation of the scalar fields using the equation (102). 
This would be regarded as a crucial difference between the AdS in light-cone and Schrodinger 
backgrounds in the finite temperature case. We defer the systematic analysis for finite M 
with full generality to the future. 

For m — 0, being required by the conformal dimension of the energy momentum tensor, 
and M = 0, being required by the condition q = (this restriction would be lifted for general 
setup), we can read off the shear viscosity from the equation (77) 

6 4 3 L s T 3 1 

" = ^sr- - ; = s- (103) 

where we use the entropy formula from [26] and take into account the factor K in front of 
the action (|9]). This is the advertised entropy- viscosity ratio of the strongly-coupled non- 
relativistic fluid advertised in [14| [16J [26]. The result is the same as the relativistic counter 
part for M — 0. 
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